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PROCESS AND APPARATUS FOR CALIBRATION OF A RADIOLOGICAL 
IMAGING DEVICE 

CROSS-REFERENCE TO RELATED APPLICATIONS 

[00011 This application claims a benefit of a priority under 35 USC 119(aHd) 

to French Patent Application No. 02 16078 filed December 18, 2002, the entire 
contents of which are hereby incorporated by reference. 

BACKGROUND OF THE INVENTION r 

[0002] This invention and embodiments thereof is directed to a medical field, 

and more particularly the field of radiological imaging devices. These devices 
comprise means for emitting radiation, such as an X ray source, and means for 
acquiring an . image, such as an image detector, mounted on the end of a C shaped 
arm, and the object of which the image is to be taken is placed between the two sides 
of the arm. The images acquired by the detector as the arm rotates around the object, 
and thus several acquisitions correspond to different points of view of the object. 
Means for processing allow a three-dimensional model of the object to be 
reconstructed firom the acquired two-dimensional images. This reconstruction 
supposes that the different positions of the device are known precisely as well as its 
geometric characteristics. . 

|0003] The three-dimensional model obtained can be used by a medical 

practitioner, such as a surgeon, before an operation in order to become familiar with 
the part of the anatomy for which the operation is intended. The three-dimensional 
model can also be used during the operation. For this, the medical practitioner 
disposes of equipment that allows the two-dimensional views of the part of the 
anatomy to be displayed in real time, these views being calculated from the three 
dimensional model. 

[0004] The reconstruction of the three-dimensional model requires the 

imaging device to be "geometrically calibrated" beforehand. This calibration allows 
the three dimensional space to be linked to the two dimensional information provided 
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by the various two dimensional projections. If the calibration is imprecise, then the 
quality of the three dimensional model reconstructed will reflect these imperfections. 

[0005] A known calibration technique comprises placing markers inside an X 

ray field, positioned on a ghost image that act as markers in space, and carrying out a 
series of acquisitions. As the position of these markers in the three dimenisional space 
is known, then the geometry of the acquisition can be deduced for each projection by 
inversion of a system of equations derived from the position of the markers on the 
projected images. A technique of this type is described, for example, in US 5 442 674. 

[0006] In general, when an image is taken of a part of the anatomy, the device 

is commanded to perform a series of acquisitions in the same geometrical conditions 
as the series of calibration acquisitions, so that the geometry of each acquisition is 
known precisely. This is possible due to the fact that the movements of the C shaped 
arm can be repeated. 

[0007] If, however, the acquisitions of the part of the anatomy are not carried 

out in the same geometrical conditions as the calibration acquisitions (number of 
views, angular positions, speed of rotation, initial and final positions), then a new 
calibration has to be performed, which increases the number of views taken for the 
calibration. 

[0008] Furthermore, the number of acquisitions required to obtain a good 

quality three-dimensional model depends on the type of anatomical structure to be 
imaged. Typically, the three dimensional reconstruction of the bone structures 
requires approximately 120 views whereas the three dimensional reconstruction of a 
blood vessel requires approximately 40 views. 

[0009] The result is that as many calibration acquisitions as possible must 

always be taken so that the calibration is valid for imaging any type of structure. 
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BRIEF DESCRIPTION OF THE INVENTION 

[0010] An embodiment of the invention is directed to minimizing the number 

of acquisitions required to calibrate the device.. An embodiment of the invention is 
directed to a process and apparatus for caHbration of a radiological imaging device 
comprising a source and a detector comprising: 

[0011] moving the device with respect to a calibration object and performing a 
series of acquisitions, each acquisition being associated to a calibration position of the 
device; _ , 

[0012] based on the acquisitions performed, determining the projection 

parameters associated to each calibration position of the device; 

[0013] for an additional position that has not been taken by the device during 

the acquisition step, determining the projection parameter values associated to this 
additional position according to the parameters associated to the calibration positions. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0014] An embodiment of the invention will be better xmderstood from the 

following description, which is given purely by way of illustration and is not 
restrictive, when read in conjunction with the appended figures in which: 

[0015] Figure 1 represents diagrammatically a C shaped arm type-imaging 

device; 

[0016] Figure 2 represents diagrammatically a C shaped arm supporting a 

source and a detector; 

[0017] Figure 3 represents the different positions taken by the source during 

the calibration acquisitions; 
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[00181 Figure 4 represents the additional positions that can be covered for the 

application of an embodiment of the process; 

[00191 Figure 5 is a block diagram schematically representing the different 

steps of an embodiment of the process; 

[00201 Figure 6 is a block diagram schematically representing the different 

steps of a variant of the embodiment of Figure 5; and 

[00211 Figure 7 is a graph illustrating the result after smoothing of the 

projection parameters calculated by an embodiment of the process. 

DETAILED DESCRIPTION OF THE INVENTION 

[00221 In figure 1, the imaging device comprises a base 1 fixed to the groimd, 

and upon which a C shaped arm 2 is mounted. The C shaped arms supports at each of 
its ends a source 3 formed by an X ray tube and an image detector 4. These elements 
are disposed so that the focus of the source 3 and the plane 6 of the detector 4 are 
diametrically opposed on the C shaped arm. The device is positioned close to a table 5 
on which a patient will be placed, the plane of the table 5 stretches between the two 
arms of the C shaped arm, which is to say between the source 3 and the detector 4. 
The C shaped arm 2 can be commanded in rotation around the X or Y axes, these two 
axes intersecting at a point O, the center of rotation of the C shaped arm. 

[0023] During the calibration phase of the device, a ghost image comprising 

markers is positioned on the table 5 between the source 3 and the detector 4, more or 
less at the level of point O. During this phase, the C shaped arm is commanded in 
rotation and the device performs a series of n acquisitions. Each image acquired is a 
projection of the markers, of the three dimensional space towards the plane 6 of the 
detector 4. 

[0024] In Figure 2, the C shaped arm 2 is shown schematically, the focal point 

S of the source 3 and the plane 6 of the detector 4 as well as an object 10 of which a 
three dimensional model is to be made. 
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[0025] The geometrical acquisition parameters of the device are, for example, 

as follows: six extrinsic parameters (three translations, three rotations) that defme the 
position in space of the C shaped arm in a tlu-ee-dimensional reference (X, Y, Z); five 
intrinsic parameters defming the projection in the detector plane: two zoom factors 
(ttu, av)along two axes of a two dimensional reference (u, v) linked to the plane, the 
co-ordinates (uo, v©) of the projection s of the focal point S of the source 3, as well as 
a deviation angle that is generally considered as nil. 

[0026] Take a point M of the object to be imaged of co-ordinates (Xm, Ym, 

Zm) in space and m its projected image of co-ordinates (Um, Vm) in the detector plane. 
This gives: 
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where R and.T are the rotation and translation matrices that define the position 
of the C shaped arm in the three-dimensional reference (X, Y, Z) and s is a scale 
factor. 

[0027] The projection matrix P is therefore a 3 x 4 matrix whose twelve 

coefficients are a result of combinations of the eleven parameters (six extrinsic 
parameters and five intrinsic parameters) and a scale factor. 

[0028] Figure 5 shows the different steps of a three-dimensional 

reconstruction process. This process comprises a calibration phase 100 beforehand 
and an acquisition phase 200 of the anatomical part to be imaged. 



[0029] In step 1 10 of the calibration phase, a ghost image is positioned on the 

table and the C shaped arm is driven in rotation. 
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[0030] Figure 3 represents the different positions xi, X2, Xn views taken by the 

source 3 around the ghost image 7 during the series of calibration acquisitions. 

[0031] In step 120 represented in Figure 5, means for processing calculate for 

each position Xi, X2, X3, ... Xn, of the source 3 a projection matrix PX (defmed by 
the relation (1) above). The projection P, associated to a position xi of the source is 
determined from the configuration of the markers on the image acquired, the position 
of the markers in space is known. 

— [0032] - —A calibration algorithm used by the. means for processing allows each 
calibration position Xi, X2, X3, ... Xn to be associated to a projection PXi, PX2, PX3, 

• • . P5Cn. 

[0033] In step 130, the means for calculation determine the geometrical 

parameters axi, ax2, ... axn of the device acquisition. For each projection Pxj, thejre are 
eleven of these parameters axi (six extrinsic parameters and five intrinsic parameters, 
as previously defined). 

[0034] In step 140, the means for processing estimate geometrical acquisition 

parameters ayi, ay2, . . . ayn-i, for these additional positions yi. 

_ [0035] As well as being shown on Figure 4, consider a series of additional 

positions yo, yi, yi, ... yn, spread over the trajectory of the source. The position yo is, 
for example, positioned before the xi position and the yn position is positioned after 
the Xn position. Each position yi(i = l...n-l) is located on the source trajectory 
between two successive calibration positions Xi and Xj+i (i<n) 

[0036] The estimation of the geometrical iacquisition parameters ayi, ay2, ... 

ayn-i, is possible given that the source trajectory is continuous. The result is that the 
variation of the acquisition parameters is also continuous. 
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[00371 The means for processing carry out an interpolation of the geometrical 

parameters axi, ax2, ... axn determined for the Xi, X2, ... Xn positions. Each geometrical 
parameter ayi is therefore a combination of the Bxi parameters: 

n 

[0038] This interpolation will only take account of a limited nimiber of 

calibration positions, less than five. For example, the means for processing can only 
take into account-the-two source. positions Xj-and.Xi+i that directly surroxmd the y\ 
position. For the yo and yn extreme positions, the means for processing perform an 
extrapolation of the geometrical parameters axi, ax2, . . . axn determined for the xj, X2,. . . 
Xn, positions. In the same way, this extrapolation will only take account of a limited 
number of calibration positions, less than five. For example, the means for processing 
can only take into account the two source positions Xj and Xj+i that are closest to the yi 
position. 

[0039] The interpolations and extrapolations carried out from two positions 

may be advantageously linear. For a greater number of positions to be taken into 
account, the interpolation functions can be polynomial, rational or another type of 
function. 

[0040] In step 150, the means for processing deduct from these parameters the 

Pyb Pya, Pyn projection matrices associated to the additional positions yo, yi, y2, 

v;,yn: 

[0041] The effect of the calibration phase 100 is to increase the number of 

calibration positions obtained with respect to the initial number n of acquisition 
positions. The imaging device is then commanded during the acquisition phase 200 to 
create an image of the anatomic structure of a patient 

[0042] In step 210, a patient is positioned on the table and the C shaped arm is 

driven in rotation. The image acquisition can be carried out for all or part of the 
calibration positions yo, xi, X2, y2, X3, ... Xn, yn of the source. 
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[0043] In Step, 220, the means for processing reconstruct a three dimensional 

model of the anatomic structures of the patient from the images acquired and the 
projections Pxi, Pyi, Px2, Py2, Px3, Pyn, Pxn, associatcd to these images. Figure 6 
shows a variant of the three-dimensional reconstruction procedure of Figure 5. 

[0044] The process is similar to that of Figure 6 except that the coefficients 
Pyo, Pyi, Pyn, of the projection matrices pyo, Pyi, Pyn, are directly determined by 
interpolation or extrapolation of the coefficients pxi, Px2, Pxn, of the pyo, Pyi, 
-pyn, matrices. --^ - - _ 

[0045] Thus at step 150, the means for processing perform an interpolation of 

the matrices pxi, Px2, Px3, . Pxn, determined for the positions xi, X2. X3 Xn at step 

120. Each Py matrix is therefore a combination of the Px matrices: 

i=l 

[0046] In the same way as above, this interpolation can only take into account 

a limited number of calibration positions, less than five. For example, the means for 
processing can only take into account the two source positions Xi and Xi+i that directly 
surround the yi position. For the y© and yn extreme positions, the means for processing 
perform an extrapolation of the geometrical parameters axi, ax2, ... axn determined for 
the xi, X2,... x„, positions. In the same way, this extrapolation will only take account 
of a limited number of calibration positions, less than five. For example, the means 
for processing can only take into account the two source positions Xi and Xj+i that are 
closest to the ys position. 

[0047] The interpolations and extrapolations carried out from two positions 

may be linear. 

[0048] The variant of the image reconstruction process shown in Figure 6 

does not require step 130 for calculating the geometrical projection parameters, as the 
interpolation is carried out directly from the projection matrices pxi, Px2, Px3, Pyn, 



8 



14XZ1 26466 



determined during calibration. Consequently, this variant is simplified in comparison 
to the process of Figure 5, which means that the processing time is shorter. 

[0049] However, the process shown in. Figure 5 permits access to the 

geometrical parameters axi, ax2, axn. In particular, this process allows these 
parameters to be smoothed out. Therefore as shown in dotted lines in Figure 5, it is 
possible to impose constraints on certain parameters taken into account by the 
calibration algorithm permitting the geometrical calibration acquisition parameters to 
be determined, 

[0050] For example, a law is imposed defining certain of the parameters taken 

into account by the algorithm. This law may comprise setting certain parameters or 
defining them by a more complex function dependent on information known 
beforehand and/or other geometrical projection parameters. For example, it can be 
imposed that the position of the projection oif the source s is always in the center of 
the detector, which is the equivalent of setting (uo, vo). It is also possible to set the 
focal distance between the focal point of the source, which is the equivalent of setting 
the zoom factors (au, av). 

[0051] The setting of one or more of these parameters can make the results 

obtained more regular and consequently eliminate the parasites they contain. 

[0052] - -..iFigure 7 illustrates the .variation of a coefficient ty of the T matrix 
defining the translation of the C shaped arm; along the y axis for a group of positions 
defined by the angle of rotation of the C shaped arm, Curve A represents the results 
obtained when no geometrical parameters are set. Curve B represents the results 
obtained when the parameters (uo, vq) have been set. 

[0053] The described reconstruction process allows the projection parameters 

to be calculated for additional source positions for which there has been no calibration 
acquisition carried out. There can, be any number of additional positions and. situated 
anywhere on the trajectory of the source depending on the requirements of the final 
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three-dimensional reconstruction. In particular, it is possible to insert as many 
additional yi positions between the Xi calibration positions as required. 

[0054] This process can be applied to any number of additional positions, 

different from the calibration positions. The effect of this is to multiply by calculation 
the number of calibration positions obtained. Consequently, this process makes it 
possible to carry out just a limited number of acquisitions. Furthermore, this process 
avoids having to cariy out new calibration acquisitions in the event of the number of 
acquisitions performed being insufficient. The projection parameters can be 
geometrical parameters characteristic of the position of the device in space (extrinsic 
parameters) or geometrical parameters characteristic of the source and the detector 
(intrinsic parameters). These parameters can also be coefficients of a global matrix 
that defines the projection of an object in the three dimensional space in the plane of 
the detector. 

[0055] One skilled in the art may make or propose various modifications to 

the structure/way and/or function and/or result for the disclosed embodiments and 
equivalents thereof without departing from the scope and extant of the invention. 
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